Doppler cooling to the Quantum limit 
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Doppler cooling on a narrow transition is limited by the noise of single scattering events. It 
shows novel features, which are in sharp contrast with cooling on a broad transition, such as a non- 
Gaussian momentum distribution, and divergence of its mean square value close to the resonance. 
We have observed those features using ID cooling on an intercombination transition in strontium, 
and compared the measurements with theoretical predictions and Monte Carlo simulations. We also 
find that for a very narrow transition, cooling can be improved using a dipole trap, where the clock 
shift is canceled. 
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PACS numbers: 37.10. De, 37.10.Gh 



Laser cooling of atoms is a technique widely used, 
often as a first cooling stage on the road to quantum 
degeneracy. In the framework of Doppler cooling, a 
moving atom is cooled because of the difference in ab- 
sorption probabilities, induced by the Doppler effect, 
between quasi-resonant red-detuned counterpropagating 
laser beams [l|, |j| . At each scattering event the momen- 
tum is, in average, modified by a recoil unit hk, with k 
the wave number. If the steady state root mean square 
(rms) of the momentum distribution. 
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is much larger than the recoil momentum. Doppler cool- 
ing can be expanded into a semi-classical theory, where 
the momentum evolution follows a damped Brownian 
motion with a pure friction force. In the low intensity 
regime the distribution is thermal-like with. 



fjp = mkbT ~ 



7mh 
80\S\ 



(4^2 



p2 



(2) 



for a standing wave in ID. S is the laser detuning, T the 
linewidth of the transition and m the atomic mass. From 
the relation ([2]), the inequality ([T]) can be reformulated 
as 
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where uj^ is the recoil frequency. Thus the semiclassical 
Doppler theory, discussed above, is valid only for atoms 
having a broad transition, i.e. fulfilling the inequality 

In the 1980s, with the progress of laser cooling and 
trapping techniques, in parallel with precise measure- 
ments of the momentum distribution, Doppler theory 
was found to be too crude to explain the steady state 



regime of atoms with complex internal structure. The 
ground state Zeeman manifold turned out to play a cru- 
cial role leading to " sub-Doppler cooling" [j-Q . However, 
Doppler theory remains valid for a two-level atom and for 
J = — > J = 1 systems. Several attempts were made to 
quantitatively test the theory of a broad transition, using 
spin polarized alkaline atoms [8|, or more recently with 
spinless, ground state, bosonic two-electrons atoms [91- 



12l | . However, most of those measurements show higher 
temperatures than predicted by the theory. Explana- 
tions suggested include spurious heating effects coming 
from spatial fluctuations of the laser intensity [lO|, and 
excited state coherences [l3|. 

Two-electron atoms remain an ideal testing ground for 
Doppler theory. In addition to broad transitions, nar- 
row intercombination transitions also exist, where the 
inequality (|3]) no longer holds. Those later cases lead 
to new features of the steady state regime, theoretically 
studied by Castin and co-authors using a full quantum 
approach of the Doppler theory [1J|. In particular, close 
to the resonance, but still on the red side of the tran- 
sition, the momentum distribution has a non-Gaussian 
shape characterized by long tails, leading to a divergence 
of the rms momentum. Hence the minimum rms mo- 
mentum -lower bounded by the recoil momentum- is red 
shifted in frequency, with respect to the 6 = — r/2 value 
for broad transitions. The use of two-electron atoms in 
laser cooling with narrow intercombination lines is grow- 
ing, and several groups already reported laser cooling on 
narrow transitions with temperature at or close to the 
recoil limit [l5|-|22| . Most of those works were carried on 
in a magneto-optical trap (MOT) where the cooling dy- 
namic and steady state regime are modified by the pres- 
ence of the magnetic field gradient [17|,ll9|. Moreover, the 
gravity plays an important role displacing the position of 
the cloud and thus changing its magnetic environment. 



In this letter, we present some experiments done using a 
dipole trap where the clock shift is canceled and we dis- 
cuss the first quantitative test of this important theory 
of quantum Doppler cooling. 

As described in [l9|, our ®^Sr cold atomic sample is 
produced as follow; after a loading and a precooling 
stage on the broad ^Sq — > ^Pi dipole-allowed transition 
(r/27r = 32 MHz, uj^/2tt = 10.6 kHz), the atoms are 
transferred to a MOT operating on the narrow ^5*0 — > '^Pi 
intercombination transition at 689 nm (r/27r — 7.5 kHz, 
uJi-/2tt = 4.7 kHz). This final cooling stage lasts for 
130 ms and leads to a cold gas containing about 2 x 10^ 
atoms at a temperature of T = 2 jiK. A few tens of mil- 
liseconds before switching off the MOT, a far-off reso- 
nant dipole trap, centered on the atomic gas, is turned 
on. This dipole trap consists of a single focused laser 
beam at 780 nm. The laser power is 120 mW for a beam 
waist of 17 /xm, corresponding to a potential depth of 
To ~ 20 fiK. The radial and longitudinal trap frequencies 
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are respectively w_l = 670 Hz and 
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8 Hz. Because of 



the weak overlap between the dipole trap and the initial 
cold cloud, at best 1% of the atoms are transferred into 
the dipole trap. 50 ms after the MOT stage, a counter- 
propagating pair of beams, aligned with the long axis 
of the dipole trap and red-detuned with respect to the 
^So — >• ^Pi transition is turned on for 450 ms (Fig. [T]). We 
carefully balance the intensity of the cooling beam mini- 
mizing the displacement of the gas center of mass in the 
dipole trap. The cloud's spatial distribution is recorded 
by absorption imaging on the broad line at 461 nm. 

A 0.3 G magnetic bias field (-B) is applied during the 
ID cooling experiment, for two important reasons. First, 
the Zeeman degeneracy of the excited state is lifted so 
that the lasers interact only with a two-level system made 
out of the 771 = — >■ 771 = transition, which is unsensi- 
tivc to residual magnetic field fluctuations. Second, the 
orientation of B, with respect to Pysoj the linear polar- 
ization of the dipole trap beam, is tuned to cancel the 
clock (or transition) shift induced by the dipole trap on 
the transition of interest [23| . The variation of the clock 
shift with respect to the angle 6 between P780 and B 
is given in Fig. ^p, and is in good agreement with the 
predicted value. The relative accuracy of the clock shift 
cancellation, for the whole trapping potential, is about 
4 kHz, i.e., below the bare transition linewidth F. At 
this precision, one can ignore the small spatial depen- 
dency of the laser detuning induced by the dipole trap 
over the confined atomic gas. The position of the reso- 
nance has been measured with a precision of ±F/2 using 
absorption spectroscopy on the cold cloud. For technical 
reasons, the linear polarization Pgsg of the cooling lasers 
is not aligned with the bias B-field. Thus the effective 
coupling intensity is reduced by a factor cos^ a ~ 1/15, 
where a is the angle between Pgsg and B. Taking into 
account this reduction factor, the effective intensity If is 
in the range 0.03 — 0.1/,, where Ig = 3/xW/cm is satu- 




FIG. 1: (a): Outline of the ID experimental set-up. The 
cold atomic gas, at the center of the picture, is held by the 
focused far-off resonance laser beam at 780 nm. The ID cool- 
ing lasers at 689 nm are superimposed on the 780 nm beam. 
Along a perpendicular axis, an imaging system on the 461nm 
broad transition records the spatial distribution of the cloud. 
A 0.3 G magnetic bias field is adjusted in angle 9 with respect 
to the polarization of the 780 nm beam in order to cancel the 
clock shift on the ID laser system, (b): Variation of the clock 
shift as a function of 6. The circles are experimental data 
points, whereas the full red curve is the predicted behavior. 
The dashed hue corresponds to the red curve shifted in 9 to 
show that the small disagreement between the experimental 
point and the predicted value is most likely due to a system- 
atic error in the angle calibration. 



ration intensity. At those intensities, the damping times 
of the momentum are long, in the range 0.02 — 0.2 s. 
As a consequence, the steady state regime is usually not 
reached in a free falling experiment, which validates the 
use of the trapping potential. 

In Fig. [2] we show the experimental mean square mo- 
mentum (black full circles) in temperature and recoil 
units. The comparison with the analytical prediction of 
the full quantum approach jl4| (green curve) shows an 
excellent agreement. Signatures of a quantum nature of 
Doppler cooling are found, e.g., the mean square mo- 
mentum has a miniinum below the recoil, and it shows 
a divergence close to the resonance, but at a detuning 
clearly larger negative than S = — F/2 predicted by the 
broad transition Doppler theory (see [Eq. ^] and the 
corresponding plot: blue, dashed curve in Fig. [5]). Spuri- 
ous heating effects reported for broad transitions [10|, |l3[ 
are not observed, most likely because we are using a two- 
level system with a long damping time. 

Several steps are necessary in order to validate the ex- 
periment and theory comparison of the preceding para- 
graph. First, the value of the mean square momentum 
can not be deduced from the spatial expansion of the 
cloud in the trap at the steady state regime, because 
of the presence of light induced collective effects. For a 
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FIG. 2; Mean square momentuni in temperature and recoil 
units as function of cooling laser detuning. The experimental 
data (black, full circles), extrapolated for A^ — >■ (noninter- 
acting gas), are compared to the full quantum approach devel- 
oped by Castin et a/. [ijj (green, solid line). The blue, dashed 
curve shows the prediction for broad transition Doppler the- 
ory. The red points (green stars) correspond to the MC sim- 
ulation performed in (without) the dipole trap. Inset: Mea- 
sured mean square momentum as function of the number of 
atoms in the trap with 5 = —20 kHz and // = 0.06/s (red 
circles) with linear fit (dashed line). 



standard 3D magneto-optical trap, the cloud would in- 
flate due to the repulsive multiple scattering force 2J| . In 
contrast, we observe a compression of the cloud because 
of the dominance of the attractive shadow force in the 
ID configuration |25l . |26[ . The mean square momentum 
is then deduced by measuring the evolution of the cloud 
over half of the trap period, after switching off the ID 
cooling laser. For that purpose, we use the relationship 
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the rms value of the cloud size in the ID harmonic trap, 
fTz{t), to the initial values (at the switching off of the 
lasers) in real space and momentum space. 

Moreover, the mean square momentum depend on the 
number of atoms, showing that collective effects also in- 
duce extra heating. Regardless of the exact origin of this 
extra heating, one example of this dependency is shown 
in the inset of Fig. [2 We extrapolate the mean square 
momentum value to the noninteracting limit (vanishing 
number of atoms) using a linear fit. The data points for 
the mean square momentum in Fig. [2]are deduced in this 
way. 

Finally, even if the trap is loose along the cooling axis 
(w|| <C Wr), it is not clear that it does not affect the cool- 
ing process. Later on we will show that the trapping 
indeed has a major impact when the transition is very 
narrow, Wr ^ F. To explore the influence of the trap, 
we use a Monte Carlo (MC) simulation comparing cool- 
ing with and without the trap. The MC simulation is 
based on a rate equation describing the scattering events 
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FIG. 3: Upper panels: Raw data, false color image of the 
atomic cloud after ID cooling in free space, and 50 ms time 
of flight. The cooling laser are along the horizontal axis. From 
left to right the laser detunings are respectively 5 = — 33 kHz, 
S — —21 kHz and 5 = — 10 kHz, and the laser intensity is 
around 0.5/s- The small horizontal asymmetry is likely due 
to imperfect balance of the cooling beams intensity, which 
were not precisely adjusted like for the cloud in the Dipole 
trap case. Middle panels: Normalized spatial distribution 
along the cooling axis extracted from the upper images. Lower 
panels: Normalized momentum distribution extracted from 
the MC simulation. The laser detuning is the same as in 
the experiment, but the simulation is performed at the low 
intensity limit. These plots show a qualitative agreement with 
the experiment at higher intensity without the added trap 
(green open circles), as well as in the trap (red dots). The 
resonance lines correspond to the vertical, dashed lines. 



where the external degrees of freedom are treated classi- 
cally. This approach, neglecting the quantum nature of 
the external degrees of freedom, is known to be consis- 
tent with the full quantum approach in free space jl4| . 
This point is also confirmed here, where the results of the 
MC simulation in free space are plotted in Fig. [2] (green 
stars). In the trap, the dynamic is more subtle because 
of the presence of the trapping force. It turns out, how- 
ever, that for the strontium intercombination transition 
with cji- ~ 0.6r the trap does not significantly modify the 
mean square momentum in the steady state (red dots in 
Fig. [2]). We will sec later that the condition Wr ^ F 
leads to different conclusions. 

Other signatures of the quantum nature of Doppler 
cooling can be found in the shape of the momentum dis- 
tribution. In the broad transition semiclassical picture, 
the momentum distribution is essentially Gaussian since 
it remains very well confined far from the two ±Sm/k 
resonance lines. With a narrow transition, a single scat- 
tering event might be enough to bring the atom out of 
resonance. As a consequence, the momentum distribu- 
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FIG. 4: MC simulation on the calcium intercombination tran- 
sition, (a): Mean square momentum in recoil units, as func- 
tion of the cooling laser detuning. The red points (green 
crosses) correspond to the MC simulation performed in (with- 
out) the dipole trap. The green, full line and the blue, dashed 
line correspond, respectively, to the analytical full quantum 
expression [ij], and broad transition limit ([Eq. ([2])]). (b): 
Normalized momentum distribution without trap (green open 
circles), and in the trap (red dots). From left to right the 
laser detunings are respectively S — — 3.5a;r, S — —2.6ijJi and 
S = — l.BciJr. Black arrows indicate those points in (a). The 
resonance conditions correspond to the dashed vertical lines. 



tion is not Gaussian anymore and shows out-of-resonance 
long tails and dips at the resonance lines IJ, l27| . A pre- 
cise measurement of the momentum distribution has been 
done for the case of free space ID cooling on clouds with 
large number of atoms, to improve the signal to noise ra- 
tio. The laser intensity was increased to 0.5/s in order to 
reach the steady state regime during the laser interaction 
time. This is a likely explanation to why there is only 
a qualitative agreement found between the experiment 
and the MC simulations done for the low intensity limit 
(see Fig. [3])- The experimental momentum distributions 
(Fig. [31 middle panels) might be decomposed into two 
domains; between and beyond the two ±Sm/k resonance 
lines (red, dashed, vertical lines). The distribution on the 
inside is in very good agreement with the MC simulation 
since it depends only on the laser detuning. For the con- 
tribution outside the resonance lines, the dependence is 
more pronounced for the experimental cases, revealing 
that extra heating is at play. 

So far, we have discussed Doppler Cooling on the nar- 
row strontium intercombination transition with F ~ w^ ~ 
kap/m. We have shown that the trap has no major im- 
pact on the steady state regime. We will now consider the 
case of a very narrow transition, i.e. T <S^ Wi- ~ kap/m, 
where the laser excitation is well localized in the momen- 
tum space. In the trap, the laser excitation frequency is 



broaden by the oscillation of the atom and occurs above 
an energy threshold corresponding to h\S\. This configu- 
ration has strong similarit y w ith the broadband cooling 
proposal discussed in Rcf. [28| . It is expected to be more 
efhcient than single frequency cooling, since an atom out- 
side the resonance lines can still be cooled and brought 
back to the central region. As a consequence the long 
tails are reduced in the trap with respect to the free space 
case. This effect is in fact evident also for the strontium 
transition (see lower- right panel on Fig. [3]). However, 
it will be more pronounced with a very narrow transi- 
tion. As an example, MC simulations were performed 
for the calcium intercombination transition (A = 657 nm, 
F/27r = 400 Hz and ajr/27r = 11 kHz), using the same 
dipole trap parameters as previously described (see Fig. 
S]) . We observed that the momentum distribution is more 
confined in the central region, leading to a reduction of 
its mean square value by more than a factor of 3, with a 
minimum closer to the resonance line. More systematic 
studies are left for future investigation. 

To conclude, we have found an excellent agreement be- 
tween cooling on the narrow intercombination transition 
in strontium, and the quantum theory of Doppler cool- 
ing developed in [IJ]. As a major feature of cooling on 
narrow transitions, the momentum distribution can be 
decomposed into a cold part at lower momenta than one 
corresponding to the laser resonance line, and a hotter 
part at higher momenta. This latter component can be 
strongly reduced by using a very narrow transition in a 
dipole trap, where the clock shift is canceled. 
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